Voronoi and Fractal Complex Networks and Their Characterization 
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Real complex networks are often characterized by spatial constraints such as the relative position 
and adjacency of nodes. The present work describes how Voronoi tessellations of the space where 
the network is embedded provide not only a natural means for relating such networks with metric 
spaces, but also a natural means for obtaining fractal complex networks. A series of comprehensive 
measurements closely related to spatial aspects of these networks is proposed, which includes the 
effective length, adjacency, as well as the fractal dimension of the network in terms of the spatial 
scales defined by successive shortest paths starting from a specific node. The potential of such 
features is illustrated with respect to the random, small-world, scale-free and fractal network models. 
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Several natural objects and phenomena can be mod- 
eled in terms of networks of basic components inten- 
sively exchanging information between themselves. Re- 
cent findings about the behavior of real networks, notice- 
ably the small- world and scale-free models[l|, motivated 
a renewal of interest in mathematical and statistical in- 
vestigations of the properties of such complex structures. 
These studies often consist of estimating fundamental 
properties of the network, such as the degree distribu- 
tion, clustering coefficient, and average length, for sev- 
eral network parameter configurations. For instance, the 
degree distribution of networks where the probability of 
new connections is proportional to the node degree has 
been verified to follow a power law [IJ , indicating scale in- 
variance, hence the name scale-free networks. While such 
investigations have led to a series of interesting findings, 
relatively little attention, with a few exceptions has 
been given to the role of spatial relationships between the 
network nodes as constrains to the network topology. In- 
deed, connections in real networks are often inherently 
related to the spatial positions of the constituent nodes 
as well as to their spatial adjacencies. Real networks 
falling into this category include but are not limited to 
telephone and electric power distribution, transportation 
systems, computer networks, spatio-temporal gene acti- 
vation in animal development, as well as the neuronal 
networks in the mammals retina and cortex. 

The following three main kinds of spatial requirements 
are considered here: (i) connect adjacent or near nodes, 

(ii) minimize the access (according to some cost such as 
time or distance) to more distant regions of the net and 

(iii) minimize the number and/or total length of edges. 
As such items are often incompatible - e.g. minimizing 
the access between any pair of nodes implies a large num- 
ber of edges, networks are often characterized by a com- 
promise between such requirements. One of the purposes 
of the current work is to investigate how the topologies 
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of complex networks are related to spatial constraints. In 
other words, as the topology of real networks is often op- 
timized to suit specific tasks, it becomes particularly im- 
portant to obtain measurements directly expressing such 
spatial/function relationships. 

Traditional random, small-world and scale free net- 
works, as well as fractal networks, are considered in this 
paper. The network nodes are assumed to be spatially 
distributed along the two-dimensional region according 
to the Poisson distribution with parameter p. The do- 
main Q, C R 2 where the network is embedded is assumed 
to be the square (0 < x < L, < y < L). The Euclidean 
distance between any two nodes i and j is represented as 
dij. Each point (x,y) of f2 is associated to the nearest 
network node, which defines a Voronoi partition of f2 Q, 
and consequently a Voronoi network. The spatial posi- 
tion of each network node i, i = 1, 2, . . . , N, is expressed 
as Pi — (xi,yi), and the area of the Voronoi cell respec- 
tively associated to this point is henceforth represented 
as Ai. Given the initial nodes, any type of network can be 
constructed by adding connections between them. The 
degree of a specific node i is ki, and the cost (or length) 
of the edge between nodes i and j, hence Xij, is given 
by the Euclidean distance between those nodes. As the 
networks considered here are not oriented, we have that 
Ajj = Xj i. Infinite length is assigned to inexistent edges 
and not taken into account for the statistics. The short- 
est path between two nodes i and j, expressed as e^.j, 
is given by the sum of the edge lengths composing that 
shortest path. 

The concept of hierarchy provides a powerful orga- 
nizational framework to model several real networks 
El 0- Such a structuring can be immediately ex- 
tended to Voronoi networks. Indeed, the tessellation of 
the network space fl provides a natural and generic basis 
for defining a broad class of spatially constrained hier- 
archical and fractal networks. The basic construction 
principle is to define a new network inside each of the 
current Voronoi cells, as illustrated in Figure ^ Each 
added hierarchy is identified by successive integer indices 
t = 0,1,..., Q, and it is henceforth assumed that the 
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FIG. 1: A Voronoi cell defined by node Pi yt -i at hierarchy t 
and a daughter network and respective Voronoi partition at 
the subsequent hierarchical level. 



nodes can only make connections inside their respective 
Voronoi cell and with the parent node, to which they are 
connected according to some statistical role, e.g. obeying 
uniform distribution with fixed parameter 7 m , as adopted 
henceforth. At least one random connection is guaran- 
teed between the parent cell and one of the nodes of the 
daughter cell. In this way, the overall network can be 
understood as a series of spatially congruent networks 
connected along the hierarchical levels. The parameter 
7 m controls the level of interaction between distinct hi- 
erarchies. In case 7 m is low, the subnetworks behave in 
an almost independent fashion. A broad variety of con- 
necting models can be adopted inside each Voronoi cell, 
including random, small- world and scale- free approaches. 
Although hybrid networks can be obtained by using dif- 
ferent connecting schemes at each hierarchical level or 
cell, the fractal networks in this article are all homoge- 
neous and random. 

Different fractal networks can be obtained by vary- 
ing the nodes distribution inside each cell and the way 
in which the connections between the nodes are estab- 
lished. It is henceforth assumed that each Voronoi cell 
along the hierarchies have the same average number of 
nodes M throughout, so that the average total number of 
nodes at level t is M t+1 . Since M is maintained through- 
out the network, cells at subsequent hierarchies become 
progressively denser, implying smaller spatial scale and 
increased level of details. The connections at any hi- 
erarchical level are uniformly established with the same 
probability 7, and it is reasonable to make 7 = j m . The 
network starts from hierarchy t = with the M nodes 
being connected with probability j m to the initial master 
node Po. As the distribution of points inside each cell at 
level t obeys a Poisson distribution with parameter pt, 
we have M = p t (A t ) and (A t ) = l/(p t -i) 0, so that 
Pt = Pt-iM = pqM 1 . Observe that A = I?. The tes- 
sellations extend up to hierarchy t = Q such as pq < 1. 
Therefore Q < —log(po — M). The average total number 
of nodes at level t is N t = ptA , so that the average total 
amount of nodes in the network is N = 1 + SE=oP*^o- 
The average node degree at any cell in hierarchy t is 



(k t ) = 7 A/. The characteristics of the network are there- 
fore completely determined by the choice of L, 7, M and 
Po- 

The self-similar nature of such fractal random net- 
works, allied to the fact that several of the properties 
inside each cell is independent of those in other cells, al- 
low convenient estimation of many network features. In 
case the distribution of a measurement w inside any cell 
at any cell in level t can be expressed in terms of t, i.e. 
Pt (w) , it is possible to estimate its distribution from level 
to b < Q by using Equation Q] 



PH = — b — E C * P *H C 1 ) 

where (q) = AopoM* is the average number of 
Voronoi cells at level t. In case Pt(w) is constant with 
respect to t and not too dispersed, as can be the case 
with (n t ) and (fc t ) under certain circumstances 0, the 
resulting distribution will lead to the same distribution 
Pt(w), now more intensely sampled. For instance, in the 
case of fractal random networks, we have that the de- 
gree distribution for N nodes will approach that of a 
random network with M nodes. In other words, the de- 
gree distribution of the fractal network does not follow 
a power law and is smaller than for a random network 
with the same overall number of nodes. The character- 
istic spatial scale at level t, hence r t , can be defined as 
the radius of the equivalent circle with area A t . Conse- 
quently, r t — (npo)~°- 5 M°- 5 ( 1 Any network property 
directly proportional to r t , as is the case with the edge 
length, shortest path and nearest neighbor distances, 
will therefore follow a power law along t with exponent 
— 0.5log(M). Such a property will also be characterized 
by unit fractal dimension with respect to St — log(r t ). 

The following measurements are suggested and used in 
the present work to characterize the spatial properties of 
complex networks in a more complete fashion: 

Edge distance: Although the average length, mean- 
ing the average shortest path between any two network 
nodes, has been traditionally used in order to charac- 
terize complex networks Q, the distribution of the Eu- 
clidean distances between two directly connected nodes 
i and j can also be provide useful information about the 
spatial characteristics of the network. The edge distance 
is henceforth expressed as Xij, with average (A). 

Effective distance: This feature applies to each pair of 
nodes i and j, defined as the sum S of the Euclidean dis- 
tances of the edges composing the shortest path between 
those nodes, in case it exists, and the Euclidean distance 
between them, i.e. E^j — S/Xij. Its maximum value of 
1 is achieved whenever the shortest path equals the Eu- 
clidean distance, implying that all edges are parallel to 
the edge between i and j. Large values of E^j indicate 
network adherence to criterion (ii). 

Adjacency: This measurement, represented as on, is 
defined for each node i as the number of direct connec- 
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tions it establishes with its n a dj spatially adjacent neigh- 
bors divided by n ad j. Thus < on < 1. Although slightly 
similar to the clustering coefficient (e.g. Q), this feature 
differs in the sense that it takes the spatial adjacency ex- 
plicitly into account. The average of this feature for all 
network nodes, i.e. (a), can be used to quantify criterion 

Network fractal dimension: While the average length 
and network diameter Q] have been traditionally applied 
to quantify the edges length distribution along the net- 
work, we present an additional measurement, related to 
the multiscale spatial dimension which is capable of 
expressing the self-similarity/fractality of the network as 
the shortest distance nodes are found while the network 
is inundated (or visited) starting from any specific node 
i. The network fractal dimension, therefore, is a func- 
tion of the spatial scale s = log (Euclidean distance), de- 
fined for each node i and represented as fi(s). First, the 
sorted sequence pi t i, . . . ,Pij, ■ ■ ■ ,Pi.N of shortest paths 
(considering the Euclidean distance and not the num- 
ber of nodes) between i and every other node j in 
the network is obtained by using traditional methods. 
The area distribution of node i can now be defined as 
hi(s) — J2j=o Ai6{s — log(pi y j)), where S is Dirac's delta 
function. The cumulative distribution of areas is given 
as &i(s) = f^ oQ hi(w)dw. In order to obtain an inter- 
polated and smoothed version of Oi(s), this function is 
convolved with a regularizing Gaussian kernel yielding 
Bi(s) = Qi{s)*g r7 {s), where g a is the normal distribution 
with standard deviation a and * stands for the convolu- 
tion operation. Now, by making Hi(s) — log(Bi(s)), the 
fractal dimension of node i as an adimensional function 
of the spatial scale is given by Equation |2 



fi{s) = dH,(s)/ds = B[{s)/B,{s) (2) 

The higher the value of /i(s), the faster the accumu- 
lated area increases along the subsequent shortest paths 
at the respective spatial scale s. Given that the areas of 
the Voronoi cells are added at once whenever their respec- 
tive distances are reached, fractal dimensions larger than 
2 can be observed. Similarly, the adoption of edge costs 
smaller than the respective Euclidean distance also leads 
to ' superfractaV dimensions, in the sense that its dynam- 
ics while it undergoes inundation exceeds the topological 
dimension of the space where the network is embedded. 
Observe that it is also possible to obtain fractal dimen- 
sions starting from a set of points instead of a single 
point. Figure [21 shows a node spatial distribution (a), the 
area distribution obtained from the respective Voronoi 
tessellation (b), the cumulative distribution (c) and the 
fractal dimension (d), all in terms of the spatial scale s. 
Such results were obtained by starting the network inun- 
dation from the leftmost point in (a). While each local 
minimum of fi(s) identifies a bottleneck along the net- 
work topology, as clearly illustrated in this example, the 
average fractal value (/) taken along all possible scales 



provides an interesting quantification of the overall spa- 
tial connectivity of the nodes in the net. 

The four diagrams in Figure [3] show the distribution of 
the above suggested measurements, taken pairwise, ex- 
perimentally obtained from several realizations of each 
network model considering L = 121 and 7 = 0.02. For 
comparison's sake, all simulations had number of edges as 
close as possible to n = 350, as well as number of nodes 
N = 140. In the case of the fractal random network, 
we used M = 5, p = 0.00035 and 7 = 0.5. For Watts- 
Strogatz, p — IO7, with the initial configuration obtained 
by directly connecting all spatially adjacent nodes. The 
Albert-Barabasi model used mo = 10 and m = 3, and 
M = 5 and po = 0.000034 were adopted for the fractal 
random networks. Distinctively segregated clusters were 
obtained in most cases, except for (d), and the fractal 
random measurements resulted generally more dispersed 
than those obtained for the other models. This is a con- 
sequence of the relatively low number of nodes M = 5 
required in order that the total number of network nodes 
results similar to those of the other simulated models. 

From figure EUa) we have that the fractal random net- 
work led to the highest overall cluster coefficients C and 
average lengths I values. The adjacencies in (b) increase 
along the following sequence: Albert-Barabasi, random, 
fractal random and Watts-Strogatz. The distinctively 
high adjacency values and low shortest distance values 
exhibited by the Watts-Strogatz model is a direct conse- 
quence of the initial adjacent configuration used in that 
case. The more intense spatial regularity of that model 
also led to high effective length values and low edge dis- 
tances in the phase space in (c) , which is also marked by 
low dispersion of the average edge distances. The frac- 
tal random network yielded effective distances slightly 
higher than those obtained for the Albert-Barabasi net- 
work, which was also characterized by small edge dis- 
tances, which reflects the larger number of small spatial 
scale Voronoi cells. Indeed, the spatially adjacent and 
spatially uniformly distributed connections in the Watts- 
Strogatz case also implied most measurements to be less 
dispersed. As shown in (d), the random and Albert- 
Barabasi networks behave similarly with respect to av- 
erage fractal dimension and node degree, whose values 
are the highest among the considered network models. 
The fractal random network led to the smallest node 
degrees and intermediate fractal dimensions. The spa- 
tially more localized connections in the Watts-Strogatz 
network, makes its inundation smoother, implying the 
lowest fractal dimension values. Except for this network, 
a positive correlation is observed in (d), reflecting the 
fact that higher node degrees tends to favor higher frac- 
tal dimensions. Another interesting point in (d) is that 
the less spatially localized connections in the random and 
Albert-Barabasi networks induced superfractal behavior. 
As a matter of fact, the relatively smaller average fractal 
dimensions presented by the Watts-Strogatz and fractal 
random models can be taken as an indication that they 
are more intensely influenced by spatial constraints, as is 
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FIG. 2: A Voronoi network (a) and its respective area dis- 
tribution (b), logarithm of cumulative areas (c) and fractal 
dimension in terms of the spatial scale s (d). 



indeed the case. 

The above results illustrate the potential of the 
suggested network measurements for quantifying and 
analysing the properties of networks as far as the possible 
influence of spatial constraints are concerned, providing 
a comprehensive integrated framework for investigating 
real complex networks. Of particular interest are the lo- 
cal/global hierarchical and topographical organization of 
mammals' cortical structures, the relationship between 
neuronal function and connectivity and shape [9j , as well 
as the spatio-temporal gene activation patterns underly- 
ing animal development. Other interesting possibilities 
are to consider fractal networks where the subsequent 
partition of the Voronoi is done with probability reflect- 
ing the degree of the parent node and/or additional spa- 
tial constraints such as varying density along Q. 
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FIG. 3: Phase space distributions of pairs of average mea- 
surements obtained from several realizations of the considered 
network models (random = o, Watts-Strogatz = o, Albert- 
Barabasi = * and fractal random = +). 



